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Abstract 

The derivation of Student's pdf from superstatistics is a mere coinci- 
dence due to the choice of the x 2 distribution for the inverse temperature 
which is actually the Maxwell distribution for the speed. The difference 
between the estimator and the variance introduces a fluctuating temper- 
ature that is generally different than the temperature of the variance. 
Fluctuations in the energy of a composite system are transferred to fluc- 
tuations in the temperature via Levy's transform. The randomization 
of the operational parameter leading to a subordinated process requires 
the same number of degrees of freedom in the directing as in the origi- 
nal process. In place of Student's pdf we find the subordinated process to 
Maxwell's speed distribution has a Pareto distribution which is interpreted 
as the density of the length of the random velocity vector. 

There is always a certain delight to see old ideas from a new perspective, but 
there is no delight to see professed 'new' ideas that are not new. The notion 
of 'superstatistics' pQ grew out of an idea that if one began with a station- 
ary Maxwell distribution and introduced a randomized operational temperature, 
then a new process could be derived which has something to do with Tsallis' 
distribution which he obtained as a stationary condition from maximizing a 
nonadditive entropy of degree-g with respect to the energy constraint The 
generalization to 'superstatistics' then consisted in replacing the exponential in 
the Maxwell distribution by e~@ H , where H is any generic hamiltonian, and 
using an arbitrary probability distribution for the conjugate (3 one could ob- 
tain a new probability distribution for H by integrating the product of the two 
probability distributions over all values of (3. We have questioned the physical 
justification of this procedure elsewhere @|. 

The idea of randomizing temperature is not new, as Beck would have us 
believe; it can be found in ref. (SJ- It is based on subordination, which can 
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be found in Feller's book UJ. If X(t) represents a brownian motion then by 
randomizing the operational time t, a variety of new processes can be derived. 
That is, if the displacements of brownian motion are governed by the probability 
density function (pdf), 

where x is the random variable and to the parameter, then exchanging their 
roles leads to a Levy pdf, 

f X0 (t) = -^e- x °/ 2t , (2) 

for the randomized time, t, where xq is the parameter. The two pdf are re- 
lated by the transform x 2 /to = x\jt. The Cauchy process is then said to be 
subordinated to brownian motion |S] 
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where the two processes are evaluated at the same time. 

The Maxwell speed, U(T), in a single dimension, and the temperature, T, 
stand in the same relation as the displacement of a brownian motion process, 
X(t), and the time, t jHj. The Cauchy distribution for the speed is therefore 
subordinated to the Maxwell distribution in one dimension [5]. 

Therefore, Beck [5] does not have to suppose that the fluctuations in inverse 
temperature, (3 follow a x 2 distribution given by his equation (2), 

i / \ n l 2 



which we would prefer to write as 
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fuM = w^[^-) (r' 2 -^-^, (4) 



r(n/2) V 2 

on the strength of equipartition, Uyfto = n. By what we shall show is an incorrect 
association of Tsallis' distribution with the subordinated process, Beck finds 
n = 2/(q — 1), where q is the degree of the nonadditive entropy, or, what 
is commonly known as the Tsallis exponent. A single mode of 10 exists at 
(3 = (1 — 2/n)Po. The single mode therefore vanishes for all q > 2. We should 
expect a qualitative change in the behavior of the process for q < 2 and q > 2. 
That no behavioral change has been predicted, already casts doubts on the 
association of half the number of the degrees of freedom with the inverse of 

(9-1). 

It follows from the Levy transform 

Pu 2 = (3 u 2 (5) 
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that the n-dimensional Maxwell speed distribution 



corresponds to the \ 2 distribution for inverse temperature, 10} , and not Beck's 
one dimensional Maxwell pdf (6)[eqn. (|14|) below]. What has happened is that, 
in a space of n dimensions, we began with a Markov process, U((3), whose sta- 
tionary transition probability is given by !©• A whole host of new processes can 
be derived by randomizing the temperature, or its inverse, [3 > 0, corresponding 
to a new random variable 2$(/3) whose pdf is l@J. The process *B(/3) is called the 
'directing' process p. 347] P- 56]. The pdf of the new process J7(5B(/3)) 
will be given by 

f°° 2 (uhinY 1 - 1 1 

f uo (u)= hWu {P)dp = - /9 , 9 w/,, / xaw -, (7) 
Jo B{n/2,n/2) (1 + {u/u Q y) n u 

where we have set (3 = 0q so as to ensure that the two processes are in thermal 
equilibium, and £?(•, •) is the beta function. Contrary to what has been claimed 
2 , is not Tsallis' canonical probability distribution [21 eqn (9)] 

/AW = / — 77: t-i/ ' V 7T- 1 + — — • (8) 
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From symmetry considerations, all odd order moments of (JBJ) are zero. In par- 
ticular, the second moment is 

which does not reflect equipartition that follows from © at all [c/., eqn l|18|l 
below ]. 

For a generic value of n, Q is a special case of the inverse beta pdf 

^ x n/2 — 1 

= B{n/2,n/2) (1 + z)"' (10) 

where a: = {u/uq)^. In economics, Hl()|l is known as the Pareto pdf, since "it 
was thought (rather naively from a modern statistical standpoint) that income 
distributions should have a tail with a density ~ Ax~ a as x — * 00, and l|lUfl 
fulfills this requirement" 6, p. 50]. 

In n = 1 dimension, J7J is the half-Cauchy pdf for a positive variate. Just 
as the Cauchy process for the displacement is subordinated to brownian motion 
in one dimension when time is randomized, the Cauchy process for the kinetic 
energy is subordinated to the Maxwellian in a single dimension when the tem- 
perature, or its inverse, is randomized 
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In n = 3 dimensions, Q is p. 95] 



8 v 2 

f n = 3 (v) = - (11) 
7r (1 + u^)' 3 



This is the density of the length of a random velocity vector v = w + + 
in 5ft 3 dimensions. It is related to the density 

f n =i\ v ) = j-,/ ! TTTi i fun! (12) 

±s(±; to — i) (1 + tr) m 

in a fixed direction, by [H] p. 32] 

/«=a(«) = "«/«=! («)> (13) 

where the prime denotes differentiation with respect to the argument. The pdf 
l|12|) has 2m — 1 degrees of freedom. For to = 1 , i|12|l reduces to the Cauchy pdf, 
while for to = 2, implying that the distribution has 3 degrees of freedom, we 
obtain l|ll|l from (|13Jl . Hence, (jll|l is the density of the length v of a random 
velocity vector in n = 3 dimensions (SI p. 95]. 

Although J2J is valid for any n > 0, Beck's 'conditional' pdf for the speed, 
(6), is our © for n = 1. For any n 1 there is an incompatibility in the 
dimensions of the two distributions. This is responsible for the seemingly close 
appearance of the pdf of the subordinated process with the Tsallis distribution. 
In other words, if Beck's (6) is meant to be the speed distribution, it does not 
reflect the dimensionality of the \ 2 pdf for /?, (0 . Rather, Beck finds that 
the subordinated process is governed by Student's t distribution, ©. The first 
factor tends to unity as n — > oo, and for every fixed u |S] 

—log 1 + ^— -> -i/3 M 2 



so that in the limit as n — > oo we have 

fp ° (u) = v1l e ~ |/W ' (14) 

This is the n — 1 Maxwcllian of which is the conditional pdf (6) that Beck 
started out with! 

What has happened is that the derivation of the Student pdf © formally 
resembles subordination, but has really nothing to do with it. Suppose that 
n + 1 random variables u and iti,«2, ■ --,Un are independent and identically 
distributed according to the normal distribution with zero mean and standard 
deviation l/y^^o- The distribution of u is given by Ijl4|l . while the distribution 

of the square root of the average of the sum of squares, v — \ — Y^—\ 
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On account of the independence of the random variables u and v, their joint 
pdf will be given by 

UMfM = r^Tf (|)"V +1)/ V-^§^W). (16) 

The probability that u/v < t is the integral of (|16J) over the region v > and 
u < to. Introducing the transformation u = xy and v — y, whose Jacobian is y, 
gives the distribution 

" rimf* (?)" /2 "»" +1,/ 7l *f *"e->"<- w . 

Its derivative yields Student's pdf 

/(t) = -4= ii-a-Ifi + L) . (it) 
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It is just a mere coincidence that the \ 2 pdf which Beck chose to represent 
fluctuations in the inverse temperature, J2J), has the same form as Maxwell's 
density for the speed \J A J27=i u i- Quite surprisingly, the Student distribution 
H17(l is independent of the variance, 1//3q. Although this is important insofar as 
testing hypotheses of the mean of a population, because it is independendent 
of the variance, is concerned, it nevertheless implies that information has been 
lost [c/., P and below]. 

The reason why Beck has /3q as a parameter in his pdf (JSJ) is that he obtained 

1 f ) / 5(«-l)/2 e -3/3(« 2 +«/A)) dB 

A L , fl,« 2 V (n+1)/2 
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using © instead of using (J2J. Had he done so, he would have obtained 

-(n+l)/2 

The two are the same if equipartition holds, 



In general, the estimator of the variance, i X]™=i u i > w * u be different from 
the variance Pq 1 ■ The estimator is /3 _1 , and this temperature will usually be 
different from the temperature of the reservoir, Pq 1 ■ In other words, the statistic 
for the Student distribution, t = wA/52i=i u ?> ^ s n0 ^ distributed as a normal 
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random variable, as it would be had the statistic been given by t = u/ ^Jn/ (3$. 
The actual variance can be completely unknown, and it suffices the sample 
variance alone in order to make statistical predictions. Nevertheless, we expect, 
that as the number of degrees of freedom increases the distribution of t will 
be very similar to that of a standard normal variable {i.e., it is a consistent 
estimator). Accordingly, we can identify the Xf in Beck's eqn (3), which are 
distributed as \ 2 i with the kinetic energies of the individual particles, and the 
left-hand side should be n/(3, and not (3 itself. 

The idea of 'superstatistics' is that the exponential pdf [7] 

is actually a conditional pdf, and specifying a normalized pdf for what was 
previously a mere parameter, say, by the x 2 pdf 

f^) = Y^fo /2 P n/2 - 1 e- 0E °. (20) 

leads to a new distribution of E which is precisely the Tsallis pdf. However, 
<|2UI) is no arbitrarily chosen pdf for the inverse temperature, but, rather the 
directing process which results when the parameter of the original distribution 
is randomized. 

Consider the usual case of a power law for the structure function jS], 

TPwn— 1 

P{E) 
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which leads to equipartition because the partition function is Z(f3) — (3 m . 
Then the process which is subordinated to the exponential pdf l|19|l has a pdf 

MD = jf MWzM * = B( ,,J, i/2)(1 figp,/^ e»> 

which is still the inverted beta pdf [c/., eqn Q]. This is because the Levy 
transform, 

[3 E = (3E , (22) 

in (|19|) produces a randomized inverse temperature which is distributed accord- 
ing to (121 )[) with n = 2m degrees of freedom. Consequently, the inverse beta 
pdf H21(l is the Fisher-Snedecor pdf with equal degrees of freedom, m = n/2, 
or what is commonly known as the Pareto distribution. The transformation to 
a new process by randomizing the operational temperature cannot change the 
number of degrees of freedom of the system. 

'F-superstatistics'^JEj is not a directing process, but, rather, a subordinated 
process that is obtained from and 10 by averaging the product of these 
densities over a common value of the kinetic energy, viz. |S] 

'*<«-/. ^*- w ^i™w 5' <23) 
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On account of the necessity that both the original and directing processes have 
the same number of degrees of freedom, (|23l) is a symmetrical inverse beta distri- 
bution, which necessarily has the same number of degrees of freedom. Comparing 
l|23|) with (|21[1 , it is easily seen that they are interchangeable on the strength of 
the Levy transformation i|22|l . Reference to a Tsallis distribution in (3 space is 
completely inappropriate. 

The symmetry of the inverse beta pdf l|23|) means that we are dealing with a 
composite system comprised of two subsystems with the same number of degrees 
of freedom. The total energy E t = Eq + E, formed from two subsystems with 
energies E and Eq, is fixed. The inverse beta pdf (|21[1 is thus transformed into 
the beta pdf, which can be written as the composition law [3J p. 80] 

rE t p2m-l 

p(Ek) = j p(E t - E)p(E) dE = 

for the structure function. 

Thus, we conclude that the statement "Of course, other distribution func- 
tions f{(3) can also be considered which may lead to other generalized statistics" 
[Til] is devoid of meaning when the original process is distributed according to 
H19(l . More importantly, as we have concluded elsewhere, "subordination can be 
considered as the probabilistic origin of power laws in physics "[31 p. 58]. 
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